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In this letter, we propose three different strategies to extract tlie weafc mixing angle (ps of the 
Bs system using penguin-mediated decays into vectors, mainly Bs — > K*''K*^, Bs — » 4>K*^ and 
Bs — > (fxj). We also provide predictions for the longitudinal branching ratio and CP-asymmetries of 
Bs — » K*^K*° using a method that combines QCD factorisation with flavour symmetries to relate 
this decay to its Bd counterpart. 

I. INTRODUCTION 

The large amount of data collected by Babar and Belle, the progress of CDF and DO and the 
advent of LHCb have increased our ability of testing the CP- and flavour- violating structure of the 
Standard Model (SM), increasing our chances of discovering New Physics. The phenomenology 
of penguin-dominated hadronic _B-dccays is particularly relevant in this fleld. At the theoretical 
level, & — > s penguin transitions are expected to receive a large impact from New Physics 
compared to 6 — *■ d penguin transitions. This comparison has driven detailed experimental 
and theoretical analyses devoted to S-decays that proceed through a. b ^ s transition, such as 
B ttK 0, B {(j),ri,uj, p . . .)Ks [1, Q ... Deviations from theoretical expectations were 
observed for some of these decays, but it remains unclear whether these results can or cannot be 
explained within the SM [i,0|- 

The properties of the Bg meson have attracted a lot of attention recently : the Bg — B^ mass 
difference AMg has been measured Q, with an immediate impact on New Physics @- Very 
recently, first experimental information on AP and the arg(— A/i2/ri2) of the Bg system has also 
been presented [Toj (see also Another important piece of information will come with the 

measurement of the mixing angle ^s- Its SM value is ^f^''^ = 2/3s = — 2A^i?f,sin7 ~ —2° [l^l and 
thus probes new CP violating phases with a high sensitivity. As far as decay modes are concerned, 
a considerable amount of theoretical work has been carried out to understand charmless Bg 
decays (for instance in Bg — > KK dec ays [Tsl . IT^ ) . The experimental program, focusing initially 
on Bg — > ttK and Bg — > KK modes [15| . will extend its scope by considering more and more 
classes of Bg decays. More observables could be studied on such details with the recent prospect 
of B or super- i? factories working at the T(55') and thus producing Bg-Bg pairs [l6j . 

On the theoretical side, non-leptonic Bg decays have been extensively studied, with a recent 
emphasis on the case of two final vector mesons (Bs — + VV), within the framework of QCD 
factorisation (QCDF) (see ref. [l^ for related work within the pQCD approach) and in 
the context of SU{3) flavour symmetries [l^ . QCDF and flavour symmetries provide different 
tools to tackle non-leptonic decays, with their advantages and shortcomings : the former is a 
systematic expansion in l/mi, but encounters difficulties with phenomenology due to power- 
suppressed hadronic effects, such as final-state interactions. The latter takes hadronic effects 
into account but may be affected by large corrections, up to 30 % for SU{3) relations. 

In ref. [l^l (see also [2l|, [l^, US); we have developed an intermediate approach that aims 
at combining SU{5) relations with QCDF-inspired input in a theoretically controlled way, in 
order to gain precision over the other approaches. In this letter we extend this approach to a 
larger class of decay modes and apply it mainly to Bg — > K*'^K*^ and partially to Bg (pK*^ 
and Bg (fxj) decays. These decays exhibit particularly alluring experimental and theoretical 
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features. For instance, on the experimental side, Bs j{*Oj{*o jg lij^ely to be measured easily 
in hadronic machines, due to its expected high branching ratio ~ 0(10"^) and its prospects 
are being analyzed in detail at present at LHCb. In addition, some decay products of the K* 
resonances are charged {7t^,K^) and are easier to identify experimentally than in the case of 
Bs K^K'^. On the theoretical side, these decays allow for an accurate extraction of the 
Bs — Bg mixing angle, with a small direct CP asymmetry in the SM that should be very sensitive 
to CP-violating New Physics. 

The aim of this letter is to provide three strategies to extract (f>s from certain non-leptonic B 
decays. It is organised as follows. In Section[ITl we consider longitudinal observables for B VV 
decays, and relate them to the observables usually obtained from an angular analysis of such 
decays. In Section IIIIl we describe a general method to extract the SM hadronic parameters 
of a i? meson decay given the branching ratio and a theoretical quantity called A, which is 
the difference between tree and penguin contributions (see [lO])- We also derive useful bounds 
for the branching ratios and direct CP asymmetries as a function of A. Section [TVl deals with 
the theoretical computation of A for B^^s K*^ K*^ , <j>K*'^ , (j>(f> decay modes : for penguin- 
dominated modes, this difference is dominated by short distances and can be computed accurately 
in QCD factorisation. In Section |Vl we exploit this theoretical information to put a bound on 
the tree pollution affecting the determination of CKM angles through the mixed CP-asymmetry. 
In Section IVTl we present general expressions to extract the CKM phases a, /3, 7 and Ps from 
hadronic penguin-dominated Bd^s decays. These general expressions can be applied to B^.s — 
K*^K*^, Bd,s (j)K*'^ and 5° — > (p4> decay modes. In Section fVIIi we use flavour symmetries 
and QCD factorisation to relate Bd j{*Oj^*o ^^^^ _^ j^-*Oj^*o observables, and exploit 
this ?7-spin symmetry to constrain the Bs — Bg mixing. Finally, we discuss the three strategies 
developed and we conclude in Section IVlIII 



II. LONGITUDINAL OBSERVABLES m B ^VV MODES 



The amplitude for a B meson decaying into 2 vector mesons can be written as 
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where Ao,+._ correspond to the amplitudes for longitudinal and transversely polarized final vector 
mesons. It is also customary to use the basis Aq || j_, where = (^4+ ± AJ)\f2. 

The vector mesons in the final state decay typically into pairs of pseudoscalar particles. A 
full angular analysis of vector-vector modes provides the following set of observables: three 
polarisation fractions /o, /_l and /|| (only two of them are independent) and their CP-conjugate 
counterparts /o,_l,|[, two phases || (again, together with a total CP-averaged branching 

ratio BR^ and a total direct CP-asymmetry ^dir- The polarisation fractions are defined as 



/o 



|AoP + |A^P-H^||| 



/o,. 



|^oP + |^±P + |^||| 



A full angular analysis is available for Bd (pK*'^ from BaBar and Belle ^ 
type of analysis is expected for Bd K*'^K*^. 



(2) 



and the same 
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We will focus in this paper on observables for the longitudinal polarisation -^d™^: 

^mbf ^^'^ -^Ar^): where only Aq occurs. These observables, free from the positive and negative 
helicity components, can be predicted with a much better accuracy. Indeed the negative-helicity 
(positive-hclicity) component of the amplitude is l/m(,-suppressed (1/m^-suppressed) because of 
the nature of the interactions involved (left-handed weak interaction, helicity-conserving strong 
interaction at high energies) [13, H^- This suppression makes longitudinal observables better 
behaved and easier to compute than transverse ones. 

Some decay channels exhibit the 1 /TO^-suppression of transverse amplitudes in a very striking 
way : the longitudinal polarisation is very close to 1, e.g. /l — 97% for B p~ ■ In such 

cases, the full observables (where Aq is replaced by the sum A = Aq + A^ +^+) coincide with the 
longitudinal ones to a high degree of accuracy. On the other hand, for penguin dominated AS* = 1 
decays, Jl can be as low as ~ 50%, so that the transverse amplitudes (or ± helicity amplitudes) 
contribute significantly to the full observables. Therefore, one must determine whether purely 
longitudinal observables can be extracted from experimental measurements. 

We start from the normalized partial decay rate of i? ^ V1V2, where the two vector mesons 
go subsequently into pairs of pscudoscalar mesons. It can be written p7j 



(PV 9 1 



Td cos Old cos 02d4> Stt |^op + + I^J 



(3) 



|Aop cos^ 6*1 cos^ 02 + l^ii p- sin^ di sin^ 6*2 cos^ ( 



+ 1 p- sin^ 6*1 sin^ 02 sin^ + Re[Ao A||] — r= sin 26*1 sin 26*2 cos 
2 2v 2 

+Im[Ao A^] — ^ sin 26*1 sin 26*2 sin (p + lm[AJ^A±]— sin^ 0i sin^ 6*2 sin 2(j) 
2\/2 2 

where {0i , 6*2 , 0) are angles introduced to describe the kinematics of the decay B V1V2 followed 
by Vi PiP{ and V2 — > ^2^2- ^1 is the angle of one of the Vi decay products in the rest frame 
of Vi relative to the motion of Vi in the rest frame of the i?- meson (same for 02 with V2). is 
the angle between the two planes formed by the decay products of Vi and V2 respectively (see 
for instance Fig.l of ref. [23| for a representation of the angles). 

There are different ways to perform the angular integrations in order to extract the purely 
longitudinal component from the differential decay rate. A first option consists in computing 
moments of cos^i (or equivalently cos(?2) : 

= / 1 a AA. f^cos^^i - dcos0^dcos02d^ ^ gps\AQ?lTB (4) 

J dcos0idcos02d(j) \2 2/ 

where gps is the product of phase-space and lifetime factors 



9 PS = ^[Ml-{m,+m2Y][Ml-{m,~m2n (5) 

A second possibility amounts to performing asymmetric integrations over one angle (28j 

plong _ 



ficos6'i / ^0036*2 / d(j)- — ^ dPsl^of /rs (6) 

_i J J- Jq dcos0idcos02d(p 
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27r/3 , 



i-sin02)d02 (7) 
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In the same way we can obtain the CP-conjugate r^°'^s{B^ /) from the corresponding CP- 
conjugate distribution, leading to the CP-averaged branching ratio of the longitudinal component 

Bi?i°ng ^ Z£ (ri°"g(i3^" ^ /) + ri°"g(s," ^ /)) = 5p5^ii^iJM (8) 

where is the CP-conjugate amplitude of ^o- 

If we include the dependence on time in the above expressions, B-B mixing modifies the 
expressions [l^. We will focus on CP-eigenstates fcp in the final state K*^K*^ and (pcj), as well 
as 0isr*° with a subsequent decay of i^r*° into a CP-eigcnstate (i^sTr" or K^tt^). 

The time evolution of these observables is obtained by considering the time dependence of 
Aq (t) [13] ■ Inserting this time dependence one arrives at the usual expression for the longitudinal 
component of the time-dependent CP-asymmetry: 

, ^ r'°"g(i3g(t) ^ /cp) - r'°"g(^°(0 ^ /cp) _ cos (AMQ + sin {AMt) 

CPl ) ri°"g(BO(t) ^ /cp) + ri°"g(BO(0 ^ /cp) cosh (Art/2) - y^^'j!,^ sinh (ArV2) ^ ^ 



where the direct and mixing-induced CP asymmetries are defined by: 

.long ^ l^oP-|^P .long ^ Im(e-^"Aglo) 



together with the asymmetry related to the width difference : 

Re(e-*'^*M*Ao) 



-^^r - |AoP-K|AoP ^^^^ 

is the mixing angle and AF = — F^. t?/ is the CP eigenvalue of the final state / (±1): 
Vk'°K'0 = Vcl><t> = 1; whereas rjx-o^ = 1 if A'*" decays into K^t^'^ and —1 if it decays into Kltt'^. 
In the latter case, the contribution from the strong process K*'^ Ktt is the same for both B 
and B decays and it cancels in the time-dependent CP-asymmetry Eq. (jlOp . which depends only 
on the amplitudes Aq and Aq. 

Finally, if the direct CP-asymmetries of all three helicity components are negligible, the lon- 
gitudinal branching ratio can be estimated very easily from: = i?i?*°*'^'/o. 



III. DETERMINING PENGUIN AND TREE CONTRIBUTIONS 

We consider a B meson decaying through b — > Dqq, with D ^ d,s, and restrict the discussion 
to the longitudinal component of the amplitude. However, the results are general and can be 
applied to any B PP, PV, VV decay (in the latter case, the relations hold for each helicity 
amplitude independently). We can parameterize the amplitudes in terms of "tree" and "penguin" 
contributions, 

Aq = Xlf>T + \[°>P , Aq = X^^T + \[^^P. (12) 

where a[^^ = V^j-iVub are combination of CKM elements, U = u,c. The penguin and tree 
contributions are defined through their associated CKM factor, and not from the topology of the 
relevant diagrams (even though in many cases, tree contributions correspond to tree diagrams). 
Such a decomposition is always possible and completely general in the Standard Model since the 
unitarity of the CKM matrix allows to recast contributions proportional to V^j^Vtb into the form 
of eq. ([T^ . We will follow the convention of calling "penguin" the piece proportional to V*jjVcb 
and "tree" the piece proportional to V*]jVub- In the particular case of penguin-mediated decays. 
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-3.15 • 10" 



-0.034 6.93 • 10" 



3.11 • 10" 



7^d 



TABLE I: Numerical values for the coefficients c[^^ and TZd for 7 = 62 



there is no actual tree diagram and the tree contribution corresponds to penguins containing a 
u-quark loop (or a t-quark loop). 

For both neutral and charged decays, one can define the CP-averaged branching ratio and the 
direct CP asymmetry as given in ([5]) and (fTUl) . From these two observables wc can obtain the 
magnitudes of the amplitudes 



We consider the quantity A defined as the difi^erence between tree and penguin hadronic 
contributions [201 



A = T- P 



(14) 



The value of A might be determined on theoretical grounds, for instance through QCD factorisa- 
tion j3l| . In the next sections, we will consider decays where such a computation is particularly 
clean and free from many long-distance uncertainties. Given the arbitrary common phase for 
T and P, we can always rotate simultaneously P and A and choose A to be real positive if we 
restrict ourselves to a given channel. We will adopt this convention in the following, unless the 
contrary is explicitly stated. 

We can write down the amplitudes p2p in the following way: 



lA, 



(D)* 



p 



P 



(15) 



The previous equations can be solved for P if we know A, and the CKM param 

eters x[f^ and xi^^ . The solutions exhibit a very simple form for A real and positive 



Rc[P] = -4^^A± 



\ 



-Im[Pl2 - 



ImfPl 



2c(^)A 



BR 

^2 



(16) 
(17) 



where the coefficients c,-^' are given by 

=Ap)|Ai^)|sin7; 4^^ = (|Ai^)p-fA(^)|A(^)|cos7)/cf); c^^^ = |A(^) + Af )p (18) 



and BR = BR^°'^^/gps. The numerical values of these coefficients are collected in Table HI Once 
P is known, Eqs. p^ . (fT7|) and (fT6|) yield the second hadronic parameter T. 

Interestingly, two consistency conditions exist between BR}°'^^, A^^^ and A, to guarantee the 
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FIG. 1: Allowed region on the B R^°''<^ - A^°.^^ plane for Ba ^ K*°K*°, according to the value of 
existence of solutions for P (the argument of the square root in Re[P] must be positive): 




with the combination of CKM factors TZd = 2|cq \/yc2 ■ The approximation for the upper 

boimd on |.4^j"^| holds up to very small corrections in the usual situation A < 0(10"^) and 
5i?long _ 0(10-6). 

The relations derived in this section apply to all charmless hadronic B meson decays, and are 
thus of quite generic nature. As an illustration, we anticipate the results of next section and 
assume that we are able to compute A accurately for the decay Bd (denoted by 

^%*K*)- Given a measured value for the longitudinal branching ratio, the quantity A^^.^* in 
Eq. ipD]) constrains the direct CP asymmetries according to Eq. fTO)) . The allowed values for the 
asymmetry are shown in Fig. [T] 

IV. THE THEORETICAL INPUT OF A = T P 

The quantity A is a hadronic, process-dependent, intrinsically non-perturbative object, and 
thus difficult to compute theoretically. Such hadronic quantities are usually extracted from data 
or computed using some factorization-based approach. In the latter case, A could suffer from 
the usual problems related to the factorisation ansatz and in particular long-distance effects. 

However, for penguin-mediated decays, T and P share the same long-distance dynamics: the 
difference comes from the {u or c) quark running in the loops Indeed, in such decays, 

A = r — P is not affected by the breakdown of factorisation that affects annihilation and 
hard-spectator contributions, and it can be computed in a well-controlled way leading to safer 
predictions and smaller uncertainties. 
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1.3 ±0.2 0.21 ±0.02 0.24 ±0.02 0.35 ±0.15 0.06 ± 0.06 0.1 ±0.2 0.218 ± 0.004 



/K-(2GeV) 




aB.^K' 








0.175 ±0.025 


0.39 ±0.06 


0.33 ±0.05 


0.221 ±0.003 0.175 ±0.025 


0.38«:J» 


0.0 ±0.3 



TABLE II: Input parameters required in QCD factorisation to compute the quantities A's and iS's 
described in the text. The masses and decay constants are given in GeV. 



rric 






/^.(2GeV) 




ai{K*) 


a2iK*) ai{K*) 




Ai^K* 37.3% 


13.2% 


0.2% 


0% 


44.2% 


0.1% 


4.6% 0.1% 


0.3% 


A'k-K' 37.5% 


12.9% 


0.2% 


0% 


44.4% 


0.1% 


4.7% 0.1% 


0.3% 





rric 


aB-^k' 

A-O 




/i.(2GeV) 


fi ai{K') a^iK") 




f^(2GeV) aaW 




Ad 


44.2% 


2.0% 






52.3% — — 




0.4% 0.7% 


0.3% 




35.0% 




0.1% 


0.7% 


58.2% 0.7% 0.1% 


0.1% 5.0% 


0% 0% 


0% 


A" 


44.1% 








52.3% — — 


— 2.1% 


0.4% 0.7% 


0.3% 



TABLE III: Relative contributions from the inputs to the errors in A for the various decays. 



For vector-vector final states, a A is associated to each helicity amplitude, but we focus on 
longitudinal quantities here. We obtain for the longitudinal A of the Bd — > K*^K*^ [Bg 
K*°K*°) decay denoted by A|-.^. (A^-.^^.): 



^K'K' I = -J^Gi \Gk' (sc) - Gk, (0)1 = (1.85 ± 0.79) x lO^^ GcV (20) 



l^K'K' I - Af'K' T^^i \^K' (sc) - Gk, (0)1 = (1.62 ± 0.69) x 10"^ GeV (21) 

where Gy = Gy — r^Gy are the usual penguin functions and Ay^y^ are the naive factorization 
factors combining decay constants and form factors (see |32| for definitions). 



A^y]>y^=^mljvX'^'''{0) (22) 



The numerical values of the used inputs are given in Table [Til The contributions to each error 
from the various sources are detailed in Table Hill For the A, as well as for the other quantities 
computed in this paper, we quote as the central value the value obtained from taking the central 
value of the inputs. To estimate the error, we vary one by one each of the inputs, compute the 
difference with the central value, then add in quadrature the resulting uncertainties. The main 
sources of uncertainties are the scale of factorisation fi, the mass of the charm quark m,,, and 
the form factor A^^^'^*. 

In a similar way, we can compute the corresponding longitudinal A for the decay modes 
Bd,s 4>K*^ and B^ (pcj): 

'K'^^C^ IG^Sc) - Gm\ = (1-02 ± 1.11) X 10-^ GeV (23) 
.,0 '^^^"^c'i|G^(sc)-G0(O)| (1.16 ±1.05) X 10-" GeV (24) 





= A 




= A 


|a;j 


= A 



't'^' Attn, 



Gi \G4,{sc) - G0(O)| = (2.06 ± 2.24) x 10"" GeV (25) 
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In the following Sections wc show how to apply the results of Sections IIIII and IIVI to the 
longitudinal contribution of penguin-dominated B — > VV modes. We will see that they can be 
used to extract the Bg — Bs mixing angle and some longitudinal observables like branching ratios 
and time-dependent CP asymmetries within the Standard Model. In particular, we outline three 
different strategics to determine the Bg — B^ mixing angle (in the SM and beyond). Indeed, 
concerning New Physics we will see that under the assumption of no significant New Physics 
affecting the amplitude, while Strategy II can detect the presence of New Physics by comparing 
the obtained cjis with (/)f ^ = 2/3^, Strategy I and III can not only detect New Physics but allow 
also for the extraction of even in the presence of New Physics in the mixing. 



V. FIRST STRATEGY TO EXTRACT 0^: BOUNDING T/P 

The b ^ s penguin-dominated decays like Bg j^*oj^*o g^^.^ principle clean modes to 
extract the mixing angle (jfs- In this section and those following, refers to the same mixing 
angle that will be measured, for instance, in the mixing induced CP asymmetry of Bg — > tpcj) 
including possible New Physics contributions in the mixing. When focusing only on SM wc will 
use the notation 0s 2/3s. 

In an expansion in powers of Xu '' /x''c \ the amplitude for the decay Bg j^*^]^*^ jg given by: 



{■■<) 



Re 



K'K* 



K'K* 



sin 7 cos (j)s -\- 



(26) 



In order to determine the accuracy of this relation, we must assess the size of the CKM-suppressed 
hadronic contribution T. Notice that this relation is valid even in presence of New Physics in 
the mixing. In the SM, one can derive from the Wolfenstein parametrisation that eq. ([H]) is of 
order (with A = Vus), and both pieces shown on the r.h.s of eg. 126]) are of this same order. 
However, despite the smallness of the ratio |Al'*-'/Ai^'' | = 0.044, a significant value of the hadronic 



ratio Re (T/P) could spoil the potentially safe extraction of sine 
in ref. [33| for B nn). The deviation from sin 0s is: 



(a similar issue was discussed 



AS{B, 



All 



A 



Re ' 



P, 



K'K' 



sm 7 cos ( 



(27) 



We want to set bounds on Re(r/P), which can be related to the inputs: 



Rein.Ro(^l±^).l + Ro(^,.l 



Re(P) A 



Re(P)2 + Im(P)2 



(28) 



Eqs. p6|) and ((T7)) show that the maximum of Rc(T/P) is reached for -4^™^ = together with 
the positive branch for Re(P) in Eq. p^ . The following boimd is obtained 



T 



Re ( - ) < 1 + I -c 



(29) 



where the lower bound for BE^°^^ and the upper bound for A must be used. In a similar way, the 
minimum of Re(T/P) occurs for Ji^^^ ~ 0, for the negative branch of Eq. for the solution 
of Re(P) 



Re ( ^ ) > 1 + ( -c(-^) - ^-{c^o^lc\''Y + {II cT) BR/A^ 



(30) 



10 20 30 40 50 

^ K*°K''°) X 10'' 



FIG. 2: Absolute bounds on AS(Bs K*°K''°) as a function of BR^°''^(Bs K*°K*°). 

where the lower bound for _Bi?'°"s ^j^^j upper bound for A must be used once again. As a 
conclusion, we obtain a range for Re(r/P) from two inputs: the branching ratio _Bi?'°"s(i?s 
K*°K*°) and Af^.^., given in Eq.(PT|). 

Using Eq. (j27|) . these upper and lower bounds on Re(r/P) are converted into a bound on 
the pollution AS{Bs K*^K*"). The latter is plotted as a hmction of the longitudinal 
BR^°"s(B, K*°K*°) in FigH 

Once a measurement of j^*oj^*o^ jg available, upper and lower bounds for (j)s are 

easily obtained. For instance, if we take as a lower bound for the branching ratio BR^°'^^{Bs 
K*°K*^) > 5 X 10-6, Fig.lHgives 0.03 < AS{Bs K*°K*°) < 0.06. In the case of a moderately 
large branching ratio x*°K*^) - (30 - 40) x 10"^, the bounds get sharper, with 

0.04 < ASiBs K*°K*°) < 0.05 and 

{A'Z^iBs ^ i^*"/?*") - 0.05) < sin0, < K°;:f (i?. ^ A'*Oa'*0) _ 0.04) (31) 

The same strategy can be applied to Bg 4)K*^ and Bg (p<f> decays 

• Take the experimental value for the longitudinal branching ratio BB}°^'^ (once available), 
and the theoretical value for A from Eqs. ([M)) or (P5|) . 

• Apply Eqs. (gn]) and to constrain the range of Re(T/P). 

• Derive the allowed range for AS* according to the equivalent of ((27)) 

• From the measured value of ^j^f^f , determine (ps from 

{A'^^ - ASraa.) < Sin 0, < {A'°^^ ~ ASrmn) (32) 

A weak mixing angle (f)s different from (f)f^^ would signal the presence of New Physics. 

Interestingly, if the longitudinal direct CP asymmetry becomes available and happens to be 
inconsistent with zero, the bounds for Re(r/P) in Eq. (|29p and ([50]) can be tightened. Eq. (pS)) 
can be exploited to derive expressions similar to Eq. and (|5D|) with a non- vanishing 
leading to stronger bounds on Re(r/P) and consequently on sin^^. 
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VI. SECOND STRATEGY : MEASURING CP ASYMMETRIES AND BRANCHING 

RATIO 



In this section, we show how we can extract mixing angles and related CKM phases in a clean 
way from experimental data, the length of two sides of the unitarity triangle and the theoretical 
quantity A. The only theoretical requirement is that the decay must allow for a safe way of 
computing A. The approach is general in the same sense as in the previous section, since it can 
be applied to any B decay into two pseudoscalars or vectors. But it yields different results for 
the four groups of decays: 



1. Bd decay through a b - 

2. Bs decay through a b - 

3. Bd decay through a 6 
eigenstate) 

4. Bs decay through a b 
eigenstate) 



d process, e.g. Bd - 
s process, e.g. Bg ~ 
> s process, e.g. Bd 



(with a subsequent decay into a CP 



d process, e.g. Bs <j)K*^ (with a subsequent decay into a CP 



As far as weak interactions are concerned, the difference between Bd and Bs decays consists in 

d and b ^ s processes differ through the CKM elements X[!^J , 



the mixing angle, whereas b 
where D = d or s. 

In the case of a Bd meson decaying through a b 
angles a and /3 from the identities: 



s processes differ through the CKM elements Ai, 

D process {D = d,s), we can extract the 



sin^ a 



BR 



2|aL^)|2|A|2 
BR 



2|A, 



(^)|2 



Al 



l-yi-(^dir) -(Anix)^ 
(^1 - V 1 - (^clir)' - (Anix)' 



(33) 

(34) 



In the case of a Bs meson decaying through a b ^ D process (D = d,s), we can extract the 
angles Ps 3J| and 7, assuming no New Physics in the decay, from the following expressions: 



sin^ Ps 



sin^ {Ps + 7) 



BR 



2IA 



BR 



2IA 



1- Vl-(Air)^-(A.ix)' 



1- Vl-(Air)'-(A.ix)' 



(35) 
(36) 



If the obtained Ps differs from its SM value, this would signal the presence of New Physics. 
Notice that this strategy is obtained by combining the definition of A with the unitarity of the 
CKM matrix, so it is designed to work only in the context of the SM. Consequently the previous 
expressions should be understood as a way of testing the SM. This is an important difference 
with Strategies I and III where one can obtain a value for the weak mixing phase also in the 
presence of New Physics in the mixing (but not in the decay). 

While the previous equations arc quite general (they can be used for B PP decays), it 
is understood that BR and Adir,mix refer to the longitudinal branching ratio and longitudinal 
CP-asymmetries, respectively, when they are applied to B VV decays. 

Eq. (|55|) provides a new way to perform a consistency test for the SM value of | sin Ps \ from 

and Si?i™s(B, ^ K*"K*"). 
using the corresponding sum 



the measurements of A^°-^^{Bs K*°K*°), A^^^{Bs ^ K*^K 
The same strategy can be applied to Bs (j)K*^ and B 
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rules. This sum rule offers several advantages : it is independent of CKM angles, and all the 
hadronic input is concentrated on a single well-controlled quantity A. 

Note that all these equations depend actually on the corresponding branching ratio and .4^p^ . 
The asymmetry is indeed related to the direct and mixing-induced CP-asymmetries through 

the equality +_L^ml?)^ + (-^Ar*^)^ 1. It was already noticed in ^ in the context of 

Bg — > K'^K~ and in [S^l in the context of i? ^ J/il)K*,D*'^D* decays that it is possible to 
extract -4^p^ directly from the "untagged" rate: 

r'™«(S,(t) ^ VV) + T^°''^(B:it) ^ VV) cx i?He-rH't + Ri^e'^^^'' (37) 
If the time dependence of both exponentials can be separated, one obtains 



A'llHBg ^ VV) = 



(38) 



The branching ratio and A^^^^ are thus the only required observables to extract /3s through this 
method, which offers the advantage of concentrating in A all the hadronic input needed to bound 
the tree-to-penguin ratio. 



VII. THIRD STRATEGY : RELATING k*°K*° AND ~* K*" K*° 

Once an angular analysis of Ed — *■ K*^K*^ is performed, it is possible to extract the CP- 
averaged branching ratio corresponding to the longitudinal helicity final state. Eqs. p6|) and (fT7|) 
can be used to extract the hadronic parameters, if one assumes that no New Physics contributes 
in an appreciable way. If flavour symmetries are sufficiently accurate for this particular process, 
this estimate can be converted into a fairly precise determination of hadronic parameters for the 
5 — > s channel Eg K*°K*°. For B^.s — *■ KK modes we noticed that J7-spin analysis 
combined with QCD factorisation led to tight constraints on the ratio of the tree contributions 
to both decay modes, as well as that for the penguins. In this section we show how to relate 
Bd x*^K*^ and Eg j(^*Oj^*o (^Q^ay modes following the same approach. 

We define the parameters 5^*^' ^"^^ ^'k'k* 

Pk'K* = fP^,K'il + 6^,^.), T^,K' = fT^,j,,a + Sl.^,) (39) 
where the factor / is given by 



0.88 ±0.19 (40) 



We compute | using QCDF. These parameters are affected by the model dependent 

treatment of annihilation and spectator-scattering contributions, so the results should be con- 
sidered as an estimate. A significant part of long-distance dynamics is common to both decays, 
and we find the following upper bounds 

|C.|<0.12, \Sl,^,\<0.15 (41) 

where the largest contribution comes from the lower value of Xb- 

We could in principle apply the same strategy to Bd,s — * (pK*^ , but the corresponding 5's are 
much larger. Indeed, the computation leads to corrections up to 5^,K^ ~ 50%. This shows that 
[/-spin symmetry cannot be expected to hold at a high accuracy for any pair of flavour-related 
processes. K^*^K'^*^ offer a much more interesting potential than other final states such as (j)K*^. 
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FIG. 3: Longitudinal branching ratio for Ba j^*Oj^*o terms of the longitudinal Bd j^*Oj^*o 
branching ratio. The light-shaded area corresponds to the uncertainty on the ratio of form factors /, 
whereas the dark-shaded area comes from varying the various hadronic inputs. 

Moreover, we cannot perform a similar analysis for (j)(f) since Bd — > is a pure weak-annihilation 
process, contrary to Bs — > (f>4> mediated through penguins. Therefore we focus on the precise 
Bs ji^*Oj^*o ]-Qodes in the remaining part of this section. Notice that the large hadronic 
uncertainties affecting B^ ifxp and Bs — > (t)K*^ have no impact when we use these modes in 
the strategies described in Sec. V and VI, since we exploited a quantity A where they cancel 
out. 

Once the hadronic parameters Pf^^j^. and T^^j^, have been obtained from Eg. ([5^ . one can 
give predictions for the Bg j^*Oj^*o observables. Note that the branching ratio BB}°^^{Bd 
K*^K*'^) is an experimental input in this analysis, and this piece of information is not available 
yet. The result for the branching ratio of Bg x*°K*'^ is given in terms of the Bd K*°K*'' 
branching ratio in Fig[3] Once the branching ratio of Bd x*'^K*^ is measured one can use 
this plot to find the SM prediction for BR^°'^^{B, ^ K*"K*"). 

The ratio of branching ratios K*^ K*°) / BE}°''^{Bd /v *°/\ *°) and the asym- 

metries turn out to be quite insensitive to the exact value of BE}°^^{Bd x*^K*^) as long as 
BR^°''s{Bd K*°K*^) > 5 X 10"'^. The numerical values are summarised in Table[IVl 

Under the standard assumption that New Physics contribution to b d penguins is negligible, 
and since the experimental input comes entirely from Bd ji[*Oj^*o (a 6 — ^- d penguin), the 
results given in Table [TVl are SM predictions. In presence of New Physics in 6 ^ s penguins the 
full prediction can be obtained by adding to the SM piece extra contributions to the amplitude 
and weak mixing angle as explained in (stI . [ssj . 

One may also use this as a strategy to extract the mixing angle 0s. If one assumes no New 
Physics in the decay Bs K*°K*°, this method relates directly A^Z^{Bs K*°K*°) and 
<Ps- FigH] shows A^^'^^{Bs j^*oj^*o^ .^g^ Once this asymmetry is measured, this plot 

can be used as a way to extract and this result can be compared to the one found in tree 
decays such as B ^ DK. A disagreement would point out New Physics. Moreover, it is 
possible to distinguish whether New Physics affects the decay or the mixing itself, by confronting 
Bi?i°"s(Bs K*°K*°) and A^^.^^{Bs K*°K*^) with the SM predictions given in Table |IV1 
If the predictions for the branching ratio and the direct CP asymmetry agree with experiment, 
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= 17 ±6 

SM 






0.000 ±0.014 




0.004 ±0.018 



TABLE IV: Results for the longitudinal observables related to Bs — > K*"K*° according to Sec. IVllI 
These are predictions for the SM contributions under the standard assumption of no New Physics in b — > d 
transition. We used <^f = 2/3, = -2° for A\°^^, and we assumed BE}°'^<^{Bi k*° K'°) > 5 x 10"''. 
The quoted uncertainty includes the errors associated to all input parameters including the variation of 
7 inside the range 56° < 7 < 68° [36[ . 




-15 -10 -5 5 10 15 

(t)s (Degree) 



FIG. 4: Relation between J^^^{Bs x*°K'^) and the Bs - Bs mixing angle 4>s. We assumed 
BR^°^^{Ba — > K*'^K*°) > 5 x 10~^ and 7 = 62°. A measurement of this asymmetry leads to a prediction 
for 0s, which includes hadronic pollution and SU(3) breaking effects, according to Sec. IVIII 

but the <f>s extracted from A^™^^{Bs — > x*'^K*'^) differs from 0f^^, this will be a clear indication 
of New Physics in Bs — Bs mixing. An interesting comparison will be allowed between the value 
for (j)s obtained here and the measurement of (j)s from the mixing induced CP-asymmetry of 
Bs ^ DK decay H^. 

VIII. DISCUSSION 

The increasing list of measured non-leptonic two-body Bd- and -B^-decays provides many tests 
of the CKM mechanism of CP violation in the Standard Model. Of particular interest is the 
determination of angles through time-dependent CP-asymmetries. For instance (pg, related to 
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Strategy 1 


Strategy 2 


Strategy 3 


Inputs 




-4d°r(S. K'°K'°) 


BB}°'^^{Bd -> K*°K*'^) 
Ajf.if,, 5t, Sp, 7 


Outputs 


(l>s 


|sin/3s|, 7 




Advantages 


Applies also to 
Bs (l)K"' and 
Bs M 


Applies also to 
Bs (l)K*° and 
^ <;/)<^ 


It can be easily 
generalized to 
include New Physics 
in the decay and 
mixing. 


Limitations 


It assumes no New 
Physics in 6 ^ s 
decay. 


It assumes no New 
Physics in 6 ^ s 
decay. 


Does not apply to 
Bs <t>K*° or 

Bs — > 4>cj> because 
5t.p are big. 



TABLE V: Comparison between the three strategies for Bs -> K'°K 



Bs — Bs mixing, should be constrained : it is tiny in the Standard Model, and can be measured 
through many penguin-dominated decays. However, for such determination to be valid, one must 
assess the size of the various hadronic quantities involved as precisely as possible. 

In the present note, we have applied ideas presented in ref. [2^ for Bd^s KK to vector- 
vector modes mediated through penguins : -B^.s — > and (j)K*^ (with the condition 
that K*^ decays into a definite-CP eigenstate). In order to combine flavour symmetries with 
QCD factorisation, we have restricted our analysis to longitudinal observables, which are under 
better theoretical control. These observables have been related to the angular analysis performed 
experimentally in Sec. |TT1 Penguin-mediated modes offer the very interesting feature that the 
difference between tree and penguin contributions A = T — P should be dominated by short- 
distance physics. It can be computed fairly accurately using QCD factorisation, and it can be 
used to determine tree and penguin contributions from observables as explained in Sees. IIIII and 
IIVI This theoretical piece of information is used to relate CP-asymmetries of Bd s — *■ K*^K*^ , 
(fxj) and (f>K*'^ to CKM angles according to different strategies. For illustration, we have focused 
on Bd,s K*'^K*°, where all three strategies apply. 

In Sec. El we have proposed to use A = T — P to put stringent bounds on the pollution due 
to hadronic uncertainties. Indeed, even though the ratio IaI^V-^c'*'' | = 0.044 is small, a large 
value of the hadronic quantity Re(r/P) could spoil the naively safe extraction of sin^g from the 
mixed asymmetry of Bg — > K*^K*'^. This strategy to control the pollution can be applied to all 
penguin-mediated processes of interest here. 

In Sec. IVI[ we have suggested a second approach, using A^^'^^{Bs K*'^K*'^), A^^^.^{Bs — > 
K*°K*°) and BP'°"s(Ss k*°K*°) to extract Isin/S^]. In principle, one can also use an alter- 
native set of experimental quantities : the branching ratio together with a direct measurement of 
the longitudinal untagged rate. The sum rule needed for the Bs is independent of the 

CKM angle 7 and the input on hadronic dynamics is limited to a single well-controlled quantity: 
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A|f,^. . This strategy can also be applied to extract f3s from Bs <f>K*'^ and Bg (fxp using 
the corresponding sum rule. 

In Sec. IVIIi we proposed a last method to determine 0s, by relying on the prediction of 
the mixing induced CP-asymmetry ^]^Sf (S^ K*°K*°) as a function of the 
A'*"A'*°) and the theoretical input A^.^.. In this strategy, tree pollution is controlled using 
the hadronic information from flavour symmetry and QCD factorisation. The outcome of our 
analysis is presented in Fig. H) This strategy requires data on Bd x*Oj^*o ^^-^ ^^ic mixing- 
induced CP-asymmetry A^^'^^{Bs K*'^K*'^). The input from Bg decay is therefore minimal : 
A^°'^^{Bs K*^K*^), while aU other inputs can be obtained from S-factories. 

A comparison among the three different strategies discussed in this paper is given in Table |Vl 
where the needed inputs are enumerated as well as the predicted observables and the range of 
validity. 

If both hadronic machines and super-B factories [39| running at T(5S') provide enough infor- 
mation on Bg-decays, it will be interesting to compare the determination from (ps following those 
methods, which rely on penguin-mediated decays, with the value obtained from tree processes 
like Bs DK [i^ . Differences between the values obtained through these two procedures would 
provide a clear hint of physics beyond the Standard Model. In such a situation, the different 
methods presented in this letter would yield very useful cross-checks for the penguin-dominated 
vector modes. 
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